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: FACILITY: MATH LIBRARY 

: ABSTRACT: 

: MTHSDSIN and MTHSDCOS are functions which return the floating point 
; sine or cosine value of their single precision floating point argu- 
; ment (radians). The call is standard call-by-reference. 

3; MTHSDSIN R? and MTHSDCOS -R7 are special routines which are the same 


as MTHSDSIN and MTH$DCOS” except a faster non-standard JSB call is 
; used with the argument in RO and ne registers are saved. 
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; MTHSDSIND and MTHS$DCOSD are functions which return the floating point 
; sine or cosine value of their single precision floating point argu- 
; ment (degrees). The call is standard call-by-reference. 
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standard JSB call is used with the argument in RO and no registers 
are saved. 
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HISTORY: 
AUTHOR: 
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; VERSION: 2 


; HISTORY: 
; AUTHOR: 
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Edit history for Version 2 


$700) - Fixed overflow problem for Large radian arguments. RNH O9-Sept=61 
-002 - pe, Uden cece for A2 = 0 in DSINCOS for small arguments. RN 
. - ept- 
6-00? - Change DSINCOS so that R6/R7 = !X! instead of X. RNH 29-Sep-81 
; 2-004 = Modified logic for converting reduced argument from integer to 
to floating format to avoid modifying the exponent of a floating 
: oint zero. RNH 21-Oct-81 
; 2-005 - Modified cosine evaluation logic to check the magnitude of YHI 
instead of YLO. RNH 01-Nov-8 
2-006 - Modified negative argument logic for DSINCOSD to eliminate bug 
; uncovered by FORTRAN QA, 
- Modified RE ye LANGE logic to fix bug detected in QAR 896. 
: RNH 14-Jan-8 
2-007 - Corrected the FFS and FFC instructions in REDUCE_LARGE to properly 
test bits 0 through 20. The loss of accuracy from only testing 20 
bits was detected in an SPR. Cleaned up some comments. JCW 8-JUN-84 
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FUNCTIONAL DESCRIPTION: 


The DSIN, DCOS and DSINCOS routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = 11*(2epi) + I*(pi/4) + YI, 


where 11 and I are integers, 0 =< I < 8 and 0 =< Y1 < pi/4. Since DSIN and 
DCOS have a period of $o0i it follows that . 


0 
DSIN(x) = DSIN(I®(pi/4) + Y1)) and 
DCOS(x) = DCOS(I*(pi/4) + Y1)). 
Ueing the trigonometric identities for the sum and difference of two angles, 
8 the following table can be generated: 
8 If l= then DSIN(x) = and DCOS(x) = 
0 DSIN(Y1) C 
1 BCOs tet! 
¢ DCOS(Y1) -DSIN( 
peimier(ee? ete 
4 -DSIN(Y1) -DCOS( 
5 pees rercenth) mei) 2 tle 
6 -DCOS(Y1) DSIN(Y1) 
7 -DSIN(pi/4-Y1) DCOS(pi/4-¥1) 
Let Y be defined as Y = Y1 if I is even and Y = if I is odd, then 


i/4 - Yi, 
each entry of the above table is of the for +/-DSINCY) or +/-DCOS(Y). Based 
on the above remarks, the DSIN, DCOS and DSINCOS routines process the input 
argument x, to obtain I and Y, and based on | selects a suitable polynomial 
approximation, p(Y), to evaluate the desired fuction. 


INPUT PARAMETERS: 


00000004 LONG =4 

00000004 x = 1*LONG ; x is input angle in radians 
00000008 sine = Gctane 3 sine is DSIN(x) 

0000000C cosine = 3*LONG 3 cosine is DCOS(x) 


OOWONA UES WN SO OD NAME WN 0 OD NAME WIN 9 OONOUE Uy 


UB RRR REESE SE FWA ROPIRIPNIPNIIIN 2 


| 
| 
| 
2: ioe Bee. ~ pcos(yt) 
4-¥1) DSINCpi /4-¥1) 
1) 
| 
| 
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| 
4; 
é ; Return sine and cosine of argument 
OOFC 8 33 ENTRY MTHSDSINCOS, “M<R2, R3, R4, RS, R6,R7> 
0 61 MTHSFLAG_ JACKET 
6D O0000000'GF 9E MOVAB G*MTHSSJACKET_HND, (FP) 
1 ; set handler address to jacket 
: 3; handler 
1 6¢ 
50 gf Bc 3070 1 6 MOVD ax(AP), RO 
00000390°EF 16 031 64 JSB MTHSDSINCOS_R7 
08 BC 0 7D 1 65 MOVQ RO, asine (AP) 
0c BC 2 7D 1 66 MOVQ R2, @cosinetAP) 
04 67 RET 
68 
69 
70 
71 -SBTTL MTHSDSIN 
48 
eR : Return sine of argument 
i 
OOFC -ENTRY MTHSDSIN, “M<R2, R3, R4, RS, Rb, R7> 


~“ 
oo 


on 
o 


MTHSFLAG_JACKET 
6D OQ0000000'GF 9E MOVAB G*MTHSSJACKET_HND, (FP) 
3 pet Senther address to jacket 
; handler 


50. 04 BC 


& 70 MOVD @x(AP), R 
00000493'EF be 


0 
JSB MTHSDSIN_R7 
RET 


-SBTTL MTHSDCOS 


; Return cosine of argument 


FUN — OOONOUS Wh 


OOF C ~-ENTRY MTHSDCOS, “M<R2, R3, R4, RS, R6, R7> 


MTHSFLAG_JACKET 


bP Ae AAD 0 tet tee te ete 


AAAA.A.NAA AA. A.A AA AAI APD PIPOPONYNUNININIPONINMINIMINIPINPONINPONONYD 


| 
, 
: 


AANA ANNI AAI IAAI IAAI AAAI AIA WIAA 


ee sss el 


SINCOS ; Floating Point Sine, Cosine ad’ pitees 16-$ AX/VMS o V04-00 pP 
ho n8d DCOS ° = SEF =}38% 9: 3: # Yate L.S REO MOS TACO.MAR: 1 - 19, 
QOO00000'GF YE is MOVAB G*MTHSSJACKET_HND, (FP) 
1 3 set handler address to jacket 
41 ; handler 
; 
70 ‘1 338 
16 4 
04 48 
034C 402 


| 
| 
MOVD ax (AP) 
400 JSB MTHSDCOS” . 
401 RET 
| 
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hh 57 00 Page ' 
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FUNCTIONAL DESCRIPTION: 


The DSIND, DCOSD and DSINCOSD routines are based on octant reduction. Given 
an argument, x, it is written in the form 


= 11*360 + 1*45 + Y1, 


where I1 and I are integers =< I < 8 and 0 =< Y1 < 45. Since DSIND and 


DCOSD have a period of 360 ie Sol tinn that 
DSIND(x) DSIND(I*45 + Y1) — 


the following table ca 


DCOSD(x) = DCOSD(1"45 + Y1) 


n be generated: 


| 

| 

eenne the trigonometric identities for the sum and difference of two angles, | 
fr) 

| 


If Il = then DSIND(x) = and DCOSD(x) = 

0 DSIND(Y1) DCOSD(Y1) 

1 DCOSD(45-Y1) DSIND(45-Y1) 

2 DCOSD(Y1) -DSIND(Y1) 

3 DSIND(45-Y1) -DCOSD(45-Y1) 

& -DSIND(Y1) -DCOSD(Y1) 

5 -DCOSD(45-Y1) -DSIND(45-Y1) 

6 -DCOSD(Y1) DSIND(Y1) 

7 -DSIND(45-Y1) DCOS(45-Y1) 
Let Y be defined as Y = Y1 if I is even and Y = 45 - Y1, if I is odd, then 
each entry of the above table is of the for +/=DSIN(Y) or +/-DCOS(Y). Based 


| 

on the above remarks, the DSIND, DCOSD and DSINCOSD routines precens the input 
argument x, to obtain I and Y, and based on I selects a suitable polynomial 
approximation, p(Y), to evaluate the desired fuction. | 
| 

| 


00000004 LONG = 4 
00000008 sind = 2*LONG 
0000000C cosd = 3*LONG 


Ft at at at at a MMMM UMMMUIUUAAUAUAUAUMUMUAUMAAtAMtHeteNatMatee 


SPE BE EWN NNN WIAD PONPIPININIPUNIND) 2 OOO OS SM DOOOOO 
MEWN OS OONAUES WN (OOD NAUE WN (OOO NOAU EWN OUOONOUS 


MT 
2 
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OOF C | ret sENTRY MTHSDSINCOSD “M<R2, R3, R4, RS, Rb, R7> 
4E 449 MTHSFLAG_JACKET 
60 00000000 ' GF 9E 8 rt MOVAB G*MTHSSJACKET_HND, (FP) 
5 3; set handler address to jacket 
0 22 ; handler 
035 
0355 450 
50 04 6C 70 0355 451 MOVD ax(AP), RO 
OO000SAB'EF 16 0359 43¢ JSB MTHSDSINCOSD_R7 
08 BC 29 7D 8 SF 4S MOVQ RQ. asind(APy 
0c BC GUS 7D 63 454 MOVQ R2, a@cosd(AP) 
0367 455 
04 0367 456 RET 
0368 457 
0368 458 
0368 459 
OOFC b3e8 rh ~ENTRY MTHSOSIND “M<R2, R3, R4, RS, R6, R7> 
036A 462 MTHSFLAG_JACKET 
6D 00000000‘ GF 9E 036A MOVAB G*MTHSSJACKET_HND, (FP) 
0371 3 set handler address to jacket 
0371 ; handler 
0371 
0371 463 
50 04 BC 70 0371 464 MOVD aX(AP), RO 
00000603'EF 16 0375 465 JSB MTHSDSIND_R7 
037B 466 
04 0378 467 RET 
037C 468 
037C 469 
037C 470 
OOFC 037C 471 -ENTRY MTHSDCOSD “M<R2, R3, R4, RS, Rb, R7> 
O37E 472 
Bere 473 MTHSFLAG_JACKET 
6D O00000000'GF 9E& O37E MOVAB G*MTHSSJACKET_HND, (FP) . 
0385 3 set handler address to jacket 
0385 : handler 
0385 
0385 474 
50 04 BC 70 OR 475 MOVD ax(AP), RO 
00000661'"EF 16 0389 476 JSB MTHSDCOSD_R7 
O38F 477 
04 O38F 478 RET 
0390 479 
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| : 1 -SBTTL MTHSDSINCOS_R7 
9 4 g ; This routine computes the DSIN and DCOS of the G-format value of RO/R1. The 
9 484 ; computation is performed one of three ways depending on the size of the 
2 : ‘ ; input argument, X: 
90 487: 1) If IX} < pi/4, then X is used directly in polynomial approximation 
3 488 ; of DSIN and D¢OS. Etat te 
90 489; 2) If pi/& =< txt < 9epi/4s, then the subroutine REDUCE_MEDIUM is called 
90 490; to reduce the argument to an equivalent argument in radians, Y, and 
90 491; the octant, I, cetera the argument. Y is then evaluated in two 
8 90 138 $ p'gnoarets chosen as a function of I, to compute DSIN(X) and DCOS(X). 
90 493; 3) If Sepi/4 =< 1X1, then the subroutine REDUCE LARGE is called to 
0390 494 ; reduce the argument to an equivalent argument in cycles, Y, and the 
0390 495 ; octant, I, contanining the argument. Y is then evaluated in two 
8395 £38 ; polynomials chosen as a function of I, to compute DSIN(X) and DCOS(X). 
0390 498 MTHSDSINCOS_R7:: 
56 50 70 th 499 move RO, R6 ; R6 = X 
OF 18 0393 500 BGEQ  POS_SINCOS ; 
QOO0039F EF 16 0395 501 JSB SINTOS : RO/R1 = DSIN(:X!), R2/R3 = DCOS(X) 
50 50 72 0398 206 MNEGD RO, RO ; RO/R1 = DSIN(X) 
05 O39E 50 RSB 
O39F 504 
039F 505 SINCOS: 
56 8000 8F AA O39F 506 BICW #*x8000, R6 3 R6/R7 = iX! 
O3A4 507 POS_SINCOS: 
56 FCS8 CF 71 Q3A4 508 CMPD D_PI_OV_4, R6 : Compare pi/4 with {Xi 
3 14 Q3A9 509 BGTR SMALC_SINCOS ; No argument reduction is necessary 
56 FCS9 CF 71 O3AB 510 CMPD 0_9_PT_OV_4, R6 : Compare 9*pi/4 with !X! 
18 0380 511 BGEQ 
00B4 = 31 bees 21§ BRW LARGE_SINCOS ; Use special logic for iXi > 9*pi/4 
0385 514; 
oye2 gl 3 pi/6 =< IKt < 9epis4 
O000069A"EF 16 0385 517 i$: JSB REDUCE_MEDIUM ; Medium argument reduction routine 
0388 518 ; R4/R7 = Y = reduced argument 
038B 8519 ; R2 = octant 
7E 54 7D Q388 520 MOVQ R4, -(SP) ; Save reduced argument on stack 
7E 56 7D OQ3BE 521 MOVQ R6, (SP) : 
52 DD Beet 5 ¢ PUSHL R2 ; Save octant bits on stack 
00000535 "EF 16 O3c3) «65 JSB M_COS ; RO/R1 = DCOS(X) 
52 E 0 bare 2 4 MOVL (SP)+, R2 ; R2 = Octart bits 
56 E D cc 5 mOVQ (SP)+, R6 : 
54 BE 7D Back 526 MOVQ (SP)+, R4 3 R4/R7 = reduced argument 
TE gpd 7D 0302 37 mova RO, - (SP) i Save DCOS(x) on stack | 
Q00004BA'EF 16 0305 528 JSB M_SIN ; KO/R1 = DSIN(X) | 
52 BE 7D DB 529 MOVQ (SP)+, R2 3; R2/R3 = DCOS(X) 
05 oF ? RSB 
Py 
O30r § 3; Logic for small arguments. ‘Xi < pi/4. | 
3 4 
3DF 5 SMALL_SINCOS: : 
56 4000 8F 861 te $ CMPW = #*X4000, RO ; Compare 1/2 with ix! : 
336 CONS E4 BLEQ 2$ ; Sufficent overghang not available 
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Cs ° : Compare with 2*-28 
7” _ a Eg 8 gerd = ™ 3 No polynomial evaluation is needed 
54 5656 § ED 540 MULD3 = R6, -R6, RO + RG/R5 = Xex 
4 Fi Fi 41 MOVQ R4, =(SP) : Put X*X on stack 
FCCE CF 5 4 73 Fa 4¢ FOLYD Ra, ACOSLENRI~1, COSTBR1; RO/RI : DcOS(x) 
ee $f 5D FD 44 MOVa RO, (SP) : Save DCOS(X) on stack 
Fo42 CF «007 5475 (0400S 545 POLYD R4, #SINLENR-1, SINTBR : RO/R1 g (x72) | 
5 88 8S Btbs a8 NDDD ORG. RO t ROURT = DSINCKS. | 
52s i it 548 MOVQ = (SP)+, R2 t R2/R3 = DCOS(X) 
oh 
70 041 551 1$: MOVD #1.0, R2 : R2/R3 = 1.0 = DCOS(X) 
50 3600 g8 AA ek 33¢ BICw #*x8600, RO + RO/R1 = 1X! 
05 0418 55 RSB 
O69 888 | 
* 9 68S 8 Otel gg CS! FIRES ORG. RO, RA Be | 
FCE2 CF «307,54 S75 0420 338 POLYD Ra, #COSLENR2-1, COSTBR2; RO/RI F acye2) 
55 FFFFIFFF ae ie a 560 BICL @#*XFFFFIFFF, RS + R4/RS = XHI 
7E 56 54 63 04 9 261 SUBD Re R6, -(SP) (SP) = XLO. mee 
- 33 ae bi er 208 MULD (SP)+, R2 > R2/R3 = XLO*(X + XHI) = A2 
08 13 0438 564 BEQL 3$ : Check for A2 = 0 
52 0080 8F A2 0430 565 SUBW #*x80, R2 > R2/R3 = A2/2 
50 52 62 0442 566 SUBD 2, RO > RO/R1 = Q(Y*2) = A2/2 
54 54 64 04452 «567 38: MULD = RG, _RG + R4/RS = XHI*2 
54 0080 BF A2 0448 568 SUBW 30s #*XB0, RG + R4/R5 = KH 22/2 
54 08 62 0440 569 SUBD 3s #1, - RG t R4/R5 = XHI*2/2 = 1 
50 540 «62045 570 SUBD 4 : RO/R1 = DCOS(X) 
52 6E) OO7D ea 571 MOVa = (SP), R2 > R2/R3 = X42 
6E 50 7D 0456 57 MOVa ~~ -RO, : Save DCOS(X) | 
FCE9 CFE «3OO7s«d52 750459 378 POLYD Re. #SINLENR-1, SINTBR ; RO/RI : 0(X72) 
36 36 $5 ace 575 ADDD _—R6,_ RO : RO/R1 = DSIN(X) 
52 BE «67D «(04665 )=— 576 MOva (SP)+, R2 + R2/R3 = DCOS(X) | 
05 0468 577 RSB 
gies a8 
0469 580 LARGE_SINCOS: ; ‘a hae | 
00000718'EF 16 Ree? 23) JSB REDUCE _LARGE : BA /aT © reduced Srpment in cycles | 
52. DD 046F 385 PUSHL R2 > Save octant bits on stack 
7E 56 D 0471 584 mova R6, -(SP) ; Save reduced 
7E 54 7D 0474 35¢ MOVQ R4, -(SP) 3 argument cn stack 
0000057C'EF 16 0477 JSB L_¢os > RO/R1 = DCOS(X) 
54 E 7D 0470 § MOVa (SP)+, R4 ; Reduced erqugent 
6 BE 7D 04 mova (SP)+, R6 : in RG/R7_ 
38 0 04 g 3 MOVL (SP)+, R2 : R2 =_octant bits 
of 8 cBs | 894 wove 80, lsh) RORY = DSTNCRD 
a BE * " F 29 mova = (SP)+, R2 t R2/R3 = DCOS(X) 
05 0492 59 RSB 
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rf 33 -SBTTL MTHSDSIN_R7 
49 a ; This routine computes the DSIN of the G-format value of RO/R1. The 
49 38 ; computation is performed one of three ways depending on the size of the 
rh 4 44 ; input argument, X: 
rh oR) 3 1) if {yi < pi/4, then X is used directly in a polynomial approximation 
3 ) . 
49 08 3 2) If pi/4 =< ix! < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
49 604 ; to reduce the argument to an equivalent argument in radians, Y, and 
49 605 ; the octant, I, containing the argument. Y is then evaluated in a 
878 606 3 pol ynoa tat chosen as a function of I to compute DSIN(X). 
49 607 ; 3) If 9epi/4 =< iKX!, then the subroutine REDUCE_LARGE is called to 
049 608 ; reduce the argument to an equivalent argument in cycles, Y, and the 
049 609 ; octant, I, contanining the argument. Y is then evaluated in a 
983 219 3 polynomial chosen as a function of I to compute DSIN(X). 
049 o1¢ MTHSDSIN_R7:: 
50 73 049 61 TSTD RO 3; Check the sign of RO 
OF 18 0495 614 BGEQ POS_SIN : 
000004A1 ‘EF 16 0497 615 JSB IN ; RO/R1 = DSIN(CiX!) 
50 50 72 Q49D 616 MNEGD RO, RO ; RO/R1 = DSIN(X) 
05 04A0 617 RSB 
04A1 618 
Q4A1 619 SIN: 
50 8000 8F AA O04A1 620 BICW #*x8000, RO : RO/R1 = {Xi 
0446 621 POS_SIN: 
50 FBS6 CF 71 Q4A6 6 § CMPD D PI_OV_4, RO : Compare pi/4 with ix! 
14 Q46AB 6 BGTR SMALC_SIN ; No argument reduction is necessary 
50 FB57 CF 71 Q4AD 624 CMPD D_9_ PT_OvV_4, RO ; Compare 9*pi/4 with iX: 
34 19 0482 625 BLSS LARGE_SIN ; Use special logic for iX! > 9*pi/4 
0484 626 
0484 627 ; 
peat ° 8 3 pi/4 =< IX! < 9epis4 
OOOO069A'EF 16 0484 630 JSB REDUCE _MEDIUM ; Medium argument reduction routine 
O4BA 631 ; R4/R7 = Y = reduced argument 
04BA O36 ; R2 = octant F 
07 OF 52 8F Q4BA 6335 M_SIN: CASEB R2, #1, #7 ; Branch to one of four polynomial 
O4BE 634 3 evaluations depending on the 
0587' 048 635 1$: «WORD P_COS_R-1$ 
0587" 400 636 -WORD 7COS"R-1$ 
06 3 4C2 637 -WORD N7SINTR-1$ 
bed * 04C4 638 «WORD N_SIN_R-1$ 
the rt) 639 «WORD N_COS_R-1$ 
0602" 04€8 640 -WORD N-COS"R-1$ 
0650" O04CA 641 WORD P_SIN_R-1$ 
065D"* O4CC ek «WORD P_SIN_R-1$ 3 octant bits. 
4CE 64 
4CE 644 ; 
O4CE 645 ; Logic for small arguments. ‘Xi < pi/4. 
4CE 666 3 
4CE 64 
4CE 648 SMALL_SIN: P 
50 3280 8F 861 #ce 649 CMPW #*Xx3280, RO ; Compare with 2*-28 | 
3 18 Q4D 650 BGEQ 1$ ; No polynomial evaluation is needed 
536 CCS 7D 0405 651 MOVQ RO, R6 3; R6/R7 = X 
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ti 
IN, 
9 8 4 0408 $34 MULD RO,RO 3; RO/R1 = Xex 
FC67 CF 5 rt$ 82 POLYD RO,  aeeaamtt SINTBR ; sd = qa(x*2) 
8 6 646 OQ4E 4 MULD R6, R 3; RO/R1 = X*q(x*2) 
6 +! GEG 655 ADDD R6, RO ; RO/R1 = DSINCX) 
4E7 $28 1$ RSB 
4E 65 
4E $28 
4E8 659 LARGE_SIN: 
00000718°EF 16 O4E 660 JSB REDUCE _LARGE ; R4/R7 = reduced argument (in cycles) 
Beee 661 ; R2 = octant bits 
54 03 4E 206 L_SIN: TSTL R4 ; Check for degenerate case 
14 1 re! o87 BEQL DEGENERATE_CASE_SIN 
07 OO 52 8 ref | 962 CASEB R2, #0, #7 
06E4' rth 667 1$ -WORD P_SIN_C-1$ 
0640° O4F 668 eWORD P_COS_C-1$ 
0640" O4FA 669 <WORD P7COS7C-1$ 
06E4" O4FC 670 «WORD P_SIN_C-1$ 
O6DA' O4FE 671 -WORD N_SIN_C-1$ 
0688' 0500 67 «WORD N_COS_C-1$ 
0688' 0502 67 “WORD N7COS~C-1$ | 
O6DA' 0504 674 «WORD N_SIN_C-1$ 
0506 675 
0506 676 
0506 677 DEGENERATE_CASE_SIN: 
0506 678 
52 8A 0506 679 BICB #1, R2 ; Compute index as (R2 - 1)/2 
52 52 FF 8F 9C 0509 680 ROTL #1, R2, R2 
03 «00 8F 050 681 SEB R2, #0, #3 
$3! ons 
0788' 0512 683 18: .WORD P_ONE-1$ | 
07CB* 0514 684 -WORD UNFL -1$ 
O7BF* 0516 685 -WORD N_ONE-1$ 
07¢B* 0518 686 -WORD UNFL -1$ 
051A 687 
| 


LT  ——— _ 
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1A 688 
1A 690 
214 24 -SBTTL MTH$DCOS_R7 
SIA 38 ; This routine computes the DCOS of the G-format value of RO/R1. The 
51A 694 ; computation is performed one of three ways depending on the size of the 
a14 $32 : input argument, X. The processing is the same as described for MTHS$DSIN_R4. 
S1A $99 : 
5073 OSA G08 tSOCOS RT RO Check f d d 
; Check for reserved operan 
50 8000 BF AA 8IC 700 BICw  #*x8000, RO  ROJRI = IX! , 
Q -FADB CF 71 1 701 CMPD DPI _Ov_4, RO ; Compare pi/4 with ix! 
1 14 05 | 1 BGTR SMALC_CO 3; No argument reduction is necessary 
50 FADC CF 71 05 70 CMPD D_9_PT_Ov_4, RO ; Compare 9*pi/4 with {X! 
47 19 3 4 Ae BLSS LARGE_Tos ; Use special logic for iX! > 9*pi/4 
O52F P86 ; 
B : ad 3 pi/& =< IK! < 9epis4 
OOOOO69A'EF 16 O52F 709 JSB REDUCE _MEDIUM ; Medium argument reduction routine 
8332 fg $ ge a : reduced argument 
$ = octan 
07 #01 52 BF 0535 ay: M_COS: CASEB R2, #1, #7 ; Branch to one of four polynomial 
0539 = 711 3 evaluations depending on the 
05D8° bets 714 1$: «WORD N_SIN_R-1$ 
0508" 0538 715 «WORD N_SIN_R-1$ 
0587" 0530 716 «WORD N_COS_R-1$ 
0587° O53F 717 WORD N~COS"R-1$ 
pees, 0541 718 «WORD P_SIN_R-1$ 
3 * 0543 719 .WORD PTSIN-R-1$ 
053C* 0545 720 «WORD P_COS_R-1$ 
053C* 0547 721 «WORD P_COS_R-1$ : octant bits. 
0549 166 
0549 as 3 
0549 724 ; Logic for small arguments. iXi < pi/4. 
0549 725; 
0549 726 
0549 727 SMALL_COS: 
50 4000 8F B61 0549 728 CMPW #*x4000, RO : Compare 1/2 with iX! 
11 14 054 729 BGTR 1$ ; Sufficent overghang is available 
56 50 7d 055 730 MOVQ RO, R6 3 R6/R7 = X 
57 FFFFIFFF 8F CA 55 731 BICL @*XFRFFIFFF, R7 3 R6/R7 = XHI 
54 50 36 $3 ssa? ; SUBD3 6, RO, R4 ; R4/R5 = XLO 
0521 1 OSSE 7 BRW NEEDS SOUBLE ; Use special logic to obtain overhang 
50 3280 8F 81 0561 734 1$ CMPW #*XxX3280, RO ; Compare with 2*- 
9A 18 266 735 BGEQ 2$ 3 No polynontal evaluation is needed 
50 0 $6 56 7 § MULD RO,RO 3: RO/R1 = X*x 
FBS7 CF 07 #850 5 298 7 POLYD RO, #COSLENR1-1, COSTBR1; RO/R1 = DCOS(X) 
05 20) 4 3 RSB 
50 08 4 378 740 2$ MOVD #1.0, RO ; RO/R1 = 1.0 = DCOS(X) 
0 57 741 RSB 
576 46 
576 74 
576 744 LARGE_COS: 
00000718"EF 16 0576 745 JSB REDUCE _LARGE ; R4/R7 = reduced argument (in cycles) 
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DWAHLMOON NE LLL SMOMOPOSLOOMOO +t 
— 
bad 


PPP PS FS OOOO 000 OO WWM ~~ 
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; R2 = octant bits 

R4 ; Check for degenerate case 
DEGENERATE_CASE_COS 


vvuse we z2v 
ANNOONM 
On OO 


Cos: 

#i, R2 
#-), R2, R2 
R2, #0, #3 
UNFL -1$ 


; Compute index as (R2 - 1)/2 


P_ONE-1$ 


MTHSOSINCOS 
2-007 


50 

OF 

50. — 8000 8F 
00000588 ' EF 
50.050 

50 =—-FA71 «CF 
000009(B ‘EF 
7E 656 

53 

00000675 'EF 
5 8E 

56 GE 

6E._ 50 
00000623°EF 
52. s«BE 


5E 10 
6— 50 
00000689° EF 
08 A 0 
E 

EF 

8E 


5 
00000637" 
52 


ao 
cov 


SSN NNN ot 
OWWOODOWOWODOOWOOONNNYYNN 2 
SOOONAVUE WN" OOONOUSWw NSD 


CGOOOCOGOOCO COCO OOOOGOCOOOCOCOCOOOCOOCOOCOOOCOOOCOOOOCOOCOOOOOOOOO 
SLU Mnnmneninoooom me Meee ee 

C009 09 G0 00 09 00 09 09 09 00 C0 00 09 C0 CO CD OD. 

NR tt 2 OOOOOCOVO 

ODOONAULS WN OVONOUS 


oN NNO NO 
MmSQoUIEoN 
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t Sine, Cosine and Sincos 16-SEP-19 
6-SEP-19 
-SBTTL MTHSDSINCOSD_R7 


The computation 
input argument, 


of DSIND and DCOSD. 


octant, 


RO 
BGEQ SINCOSD 
BICW #*x8000, RO 
JSB SINCOSD 


eee 


D.45, RO 
BGTR SMALL SINCOSD 
JSB REDUCE _DEGREES 


mova RG, -(SP) 
PUSHL R 

JSB EVAL_COSD 
MOVL (SP)F, R3 
MOVQ (SP), R6 

MOVQ RO, (SP) 

JSB EVAL_SIND 
MOVa (SP)#, R2 
RSB 
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This routine congutes the DSIND and DCOSD of the D-format value of RO/R1. 
7 performed one of two ways depending on the size of the 


1) If ixi < 45, then X is used directly in polynomial approximation 


2) If 45 =< txt, then the subroutine REDUCE_DEGREES is called to reduce 
the argument to an equivalent argument in degrees, Y, and the 

containing the argument. Y is then evaluated in two 

polynomials chosen as a function of I, to compute DSIND(X) and 


Compare 45 to iX! 

special processing for small arg 
R6/R7 = reduced argument 

R3 = octant 

Save reduced arg 

Save octant bits 

RO/R1 = DCOSD(Y) 


R6/R7 = reduced argument 
Save DCOSD(Y) 

RO/R1 = DSIND(Y) 

R2/R3 = DCOSD(Y) 


Allocate 4 longwords on stack 
Save argument 

RO/R1 = DCOSD(:X!) 

Save DCOSD(iX:) 

RO/R1 


/ = argument 
RO/R1 = DSIND(X) 
R2/R3 = DCOSD(:X!) 
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° § -SBTTL MTHSDSIND_R7 
6 4 ; This routine computes the DSIND of the D-format value of RO/R1. The | 
6 5 ; computation is performed one of two ways depending on the size of the input 
2 § 3 argument, X: | 
° 8 ; 1) If Kt 45, then X is used directly in polynomial approximation | 
6 os 2) If 45 =< txt, then the subroutine REDUCE_DEGREES is called to reduce 
6 13 the argument to an equivalent argument in degrees, Y, and the 
6 § ; octant, I, containing the argument. Y is then evaluated in two 
? z 3 polynomials chosen as a function of I, to compute DSIND(X). 
6 5 MTHSDSIND_R7:: | 
50 6 § TSTD RO 3; RO/R1 = X 
OF 06 3 BGEQ POS_SIN : 
we” 06 3 JSB NEG_SIN 3 RO/R1 = DSIND(iX!) 
50 0 06 83 MNEGD RO, RO 3 RO/R1 = -DSIND(iX!) 
be re RSB 
3 bee NEG_SIND: | 
50 8000 8F 06 se BICW #*x8000, RO 3; RO/R1 = ix! 
06 44 POS_SIND: 
50 FA16 CF 06 845 CMPD D.45, RO ; Compare 45 to iX! 
1 06 846 BGTR SMALL_SIND ; special processing for small arg 
000009(B' EF 06 847 JSB REDUCE_DEGREES ; R6/R7 = reduced argument 
06 848 ; R3 = octant | 
be 849 
6 850 EVAL_SIND: 
07 OO 53 A 06 g3! CASEB R3, #0, #7 
068E' 06 3 1$ .WORD P_SIN_D-1$ 
05E1° 0629 85 «WORD P_COS_D-1$ 
05E1' 0628 854 -WORD P7COS"D-1$ 
068E* 062D 855 «WORD P_SIN_D-1$ 
068B' 062F 856 «WORD N_SIN_D-1$ 
ee 0631 857 -WORD N_COS_D-1$ 
O62E 0633 858 "WORD NZCOS"D-1$ | 
0688" 0635 59 «WORD N_SIN_D-1$ 
637 360 
637 61 
637 6¢ SMALL_SIND: 
50 FAOS CF 71 0637 6 CMPD D_SMALLD, RO : Compare 180/pi*2*-27 with ix: 
14 tt 64 BGTR 1$ ; No polynomial evaluation is 
56 D 88 ‘ 65 MOVa RO, R6 3 )=o necessary 
0671 31 64 66 BRW P_SIN_D ; 
50 ; 44 67 1$ TSTD RO : Check for zero 
18 = «1 ret 68 BEQL $ ; Return if RO = 0 
50 FAI4G CF 71 064 69 CMPD B SMALLEST_DEG, RO ; Check for possible underflow on 
93 15 0640 870 BLEQ $ 3 conversion to radians 
06 31 064F {a BRW UNFL 3: Underflow will occur on conversion 
52 50 FOFA CF 65 o3¢ fs 2$ MULD3 D_CONVERT, RO, R2 3: R2/R3 = (pi/180 - 2*-6)*ix! 
50 0300 8F A2 065 7 SuBW «= #™X300, RO > RO/RI = 1Ki#24= | 
50 2 g 65D 74 ADDD R2, RO 3 RO/R1 = DSIND(IX!) = (pi/180)ix: 
660 75 3$ RSB 
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$ } a -SBTTL MTHSDCOSD_R7 
1 8 ; This routine computes the DCOSD of the D-format value of RO. The computation 
1 0 ; is performed one of two ways depending on the size of the input argument, xX 
0d 1 ; Details are given in the discussion on MTH$DCOSD_R4. 
$6] é MTHSDCOSD_R7:: 
30 73 066 4 TSTD RO ; Check for reserved operand 
28 8000 8F AA 066 5 BICW 4 pee RO 3 RO/R1 = iX! 
F9C4 CF 071 (066 $ CMPD 45, RO : Compare 45 to {Xx} 
A 14 066D BGTR SMALL_COSD : 
OOOOO9CB*EF 16 oof $8 JSB REDUCE_DEGREES ; R6/R7 = reduced argument 
of 4 ; R3 = octant 
67 91 EVAL_COSD: 
-T &@ FB F of 35 CASEB R3, #0, #7 
05 F* 0679 893 18: .WORD P_COS_6-1$ 
6 ie 678 94 «WORD P_SIN_D-1$ 
be * 0670 95 «WORD N_SIN_D-1$ 
OSDD' 067F 896 -WORD N7COS"D-1$ 
O5DD° 0681 9 «WORD N_COS_D-1$ 
0639" 0683 98 -WORD N_SIN_D-1$ 
38 C* 0685 99 «WORD P_SIN_D-1$ 
058F* 0687 900 -WORD P7COS"D-1$ 
83 a 
O69 308 SMALL_COSD: ; : 
50 F9B3 CF 71 0689 904 CMP D_SMALLD, RO ; Compare 180/pi*2*-27 with ix! 
14 068E 905 BGTR 1$ 3; Check if polyinomial evaluation is 
56 50 7D 0690 906 MOVQ RO, R6 3 necessary 
0572 31 0693 907 BRW P_¢oS_D : POLY needed 
50 08 ny 0696 ane 1$: nove #T, RO ; RO = 1. = DCOSD(IX!) 
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-SBTTL REDUCE_MEDIUM 


; This routine assumes that the absolute value of the argument, X, is in RO/R1 
; and that pi/4 =< IX! < 9epi/4._ It returns a pair of d-format values for the 
: ene a aes YHI in R6/R7, and YLO in R4/R5S. The octant bits in are 

3 returned in R2. 


1 
313 
Hf 
12 
it 
313 
920 
69A 921 
69A 9 § The reduced ergunent is obtained by locating the octant that X is in through 
ooh 4 Z ; @ binary search and then subtracting off a suitable multiple of pi/2 
en 3 ; REDUCE_MEDIUM 
50 8000 8F AA 069A 9 $ ~ BICW #*x8000, RO 3 RO/R1 = {Xi 
50 F975 CF 71 O69F 928 CMPD D.5_PI_OV_4, RO 3 
11 15 06A4 859 BLEQ 3 Mt >= Sepi/s 
50 F966 CF 71 Q6A6 930 CMPD 4 3_PI_Ov_4, RO 
05 15 O6AB 931 BLEQ $ 3 IX! >= 3episs 
52 01 DO O6AD 236 MOVL #1, R2 ; First quadrant 
16 «11 0680 933 BRB SUBTRACT 
52 03 00 O62 935 3$: MOVL #3, R2 3; Second quadrant 
11. 11 «0685 © 936 BRB SUBTRACT 
0687 937 
50 F965 CF 71 0687 938 5$: CMPD D_7_PI_OV_4, RO ; 
0S 15 O068C 939 BLEQ 7$ 3 IML >= 7epi/4 
52 05 0d0 O6BE 940 MOVL #5, R2 3; Third quadrant 
05 «(11 6c] 4) BRB SUBTRACT 
52 07 DO O06C3 943 7$: MOVL #7, R2 3; Fourth quadrant 
00 11 0666 944 BRB SUBTRACT 
06C 945 
06C 946 
06C 947 SUBTRACT: 
53 FD A24 3E 06C 948 MOV oo cre) C Re) R3 ; R3 = index into PI_OV-2 table 
53 F996 CF4 of bess 949 MOVA PI_OV_2CR 4, R3 ; R3 =_pointer into PI_Ov_2 table 
50 83 63 06D 950 SUBD (R3)+> RO, 3 R6/R7 = 1st approximation to YHI 
96 15 06D7 951 BLEQ 1$ Peres ee 
5 D6 0609 826 INCL R2 ; Adjust octant bits 
54 56 8000 BF AB 06DB 953 1$ BICW3 #*X8000, R6, R4 ; R4 = high 16 bits of {YHI"} 
54 2700 8F Bi O06E1 954 CMPW #*x2700, ; Check for at least 6 significant bits 
0c 16 GES 935 BGTR | NOT_ENOUGH_BITS : 
54 56 7D O6E8 957 mMOVQ R6, R4 ¢ R4/RS = YHI' 
57 D& OQO6EB 958 CLRL R7 3 R6/R7 = high 24 bits of YHI" = YHI 
54 56 66 6ED 959 SUBD R6, R4 3 R4/RS = Low bits of YHI' 
54 636 FO 960 SUBD (R$), R4 > R4/R5 = YLO 
0 6F 961 RSB 
6F4 306 
6F4 963 NOT_ENOUGH BITS: 
50 63 7D Q6F4 964 mova (R3),_RO : 
50 003F0000 8F CA O6F7 965 BICL  #*x063F0000, RO ; | 
51 D4 ore 966 CLRL R1 3 
54 83 50 =663 0 967 SUBD3 RO, (R3)+, R4 3 
6 50 62 0704 968 SUBD RO, R6 3 R6/R7 = YHI 
| 


53 D 8F 
FFFFFFES 8F CA 
QOOOO000'EF ODE 


. Ff 
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00004000 
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This routine is used to reduce large arguments (iX! >= 9*pi/4) modulo pi/4. 
It returns the reduced argument, Y, in R4/R7 in units of cycles, and returns 
the octant bits, I, in R2. 


The method of reduction is as follows: 


x*(4/pi) gape recy aeae where n is an integer and 1/2 =< f < 1 
*(m=56) #(2*56%f) #(4/pi) 

(2*°56%f) #(2* (n=56) #4/pi) 

= K*C, where K = 2°56*f is an integer and C = 2*(n=-56)*4/pi 
Let L = K*C modulo 8, where 0 =< L < 8, and let I = the integer(L) and 
h = fract(L), then if I is even Y = h, otherwise Y = 1h 


CONSTANTS: 


weights exponent by 61 

weights oxpynens by 32 

maximum unbiased exponent 

Used to locate hee point in 
MTHSAL_4_OV_PI table 

D_2_T0_32: 


QUAD *x5080 3 2°32 


~ 
m 
] 
= 
= 
m 
ol 
a 
— 
nhunnu 
> 
Oo 
ooow 
ooo 


REDUCE _LARGE: 


The first step is to obtain the location of the binary point in the represen- 
tation of C = 2“(n-56)*(4/pi) in two parts - the number of longwords from 

the start and the number of bits from the most significant bit of the next 
longword. Also K = 2°56*f must be obtained. 


BICW #*x8000, RO RO/R1 = {Xi 

ROTL #-7, RO, R3 Shift exponext field 7 bits right 

SUBW = #W_ADJUST, R3 Unbias exp and adjust for leading 

zeroes. R3 = location of binary 


point 

Divide R3 by 32 and mull by 4 to get 
R4 = # of longwords (in bytes) to 
binary point. 

ms pees address of MTHS$AL_4_OV_PI 
able 

Re geznts to Ist quadword of interest 

R3(7:0) = # of bits within longword 


#-3, R3, R4 
BICL #*XFFFFFFE3, R4& 
MOVAL MTHSAL_4_OV_PI, R2 


R4, R2 
BICB #*XEO, R3 


BICW #*X7F80, RO Clear exponent field 
W RO = 2*24ef 


BIS #*x4C00, RO : a : 

CVTDL = RO, RO ; RO = High 24 bits of K 

ROTL #16, R1, RI ; Rl = Low 32 bits of kK 

BGEQ 1$ ; Check for high bit of R1 set 
INCL RO ; Adjust RO if R1 is negative 


: The next step is to generate an approximation to C, call it C'* to be usec 
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; in comput ing weet" C'' will consist of the first three integer bits of 
; C and the first § fraction bits of C. These bits will be obtained from a 
; constant stored in the interger array MTH$AL_4_OV_PI. 


; NOTE: The ASHQ, ADDL, and MULL instructions in the follow sections may 
3 result in an integer overflow trap. The overflow incurred is intentional, 
j 


hep 


; SO that the IV bit must be turned off. The IV bit is not restored until 
; after all of the necessary fraction bits hav been generated. 


$: MOVPSL =-(SP) ; Put current PSL on stack 
BICL #*C<PSLSM_IV>, (SP) ; (SP) = current IV bit 
BICPSW #PSLSM_IV ; Clear integer overflow bit 


; The necessary calculation to produce the reduced argument can require up to 
; nine longwords of temporary work space. This work specs will be allocated 
; on the stack. The work space will be accessed though the use of three 

; registers: R4, R5, and SP. For the purposes of comments the temporary work 
; space will be referred to as locations [0 though T8. The stack and its 

; pointers will look something Like this: 


7E OC 
6E FFFFFFDF 8F CA 
2089 


TO: ! }:' een. 
ie To maaan 

Td: | ' 

T3: mere Ne tielasint 
T4: } 
TS: | 
Té6: 
T?: 
T8: 


<---= RG 


<---- RS 


The following code allocates the storage and sets up the pointers. 


#36, SP : Allocate 9 Longuords on the stack 
ADDL #16, SP, R4 : R4 points to T 
ADOL #32, SP, R5 ; RS points to T8 


Get C** = C€(O0):C(1):C(2):C(3) im T5/T8. (C(O) though (C(3) are unsigned | 
Lntegers generated from the etnary representation of C. The high three bits 
of C(O) are the the first three bits to the left of the binary point of C. 
The remaining bits C(0) and C(1) though C(3) are the first 125 bits to the 
right of the binary point of C. Note that the C(i)'s are adjusted to 
compensate for their signed (rather than unsigned) interpretation in the EMUL 
instruction. Note also that the representation of C has no more than 15 
consequtive ones, so that no carry is possible from the adjustment. 


dd ed od dd 8 dd 8 = 3 2 = 


SO 


D 14 
MTHSDSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP=-1 01: AX/VMS Macro V04-00 Page 26 | 
O~007 REDUCE _LARGE 6-SEP-19 38 § 11: £33: Hy EMTHRTL. St SMTHDSINCO.MAR: 1 (19) | 
076A 1089 ; | 
76A 1090 
7, wm OO ten ! 4 ADDL3) #12, R4, R7 ; inital tae loop counter. R7 points 
or gw b78 Hi ASHQ R3, (R2), (R7) ; shift, the proper quadword so that 
a ys O77 1095 SUBL #4, R7 ; Ry points to T6 
on ae as 1938 2$ SUBL #4, R2 ; ne poross touch quanear’ in 
rT &t O78 1998 ASHQ R3, (R2), (R7) : at 1) quadword” so that C(n) is in 
: :.8 
03 18 O77C 1100 BGEQ $ i Check tor high bit’ cof C(n) set | 
08 A? 06 Ee: 1193 INCL ate?) ; Bit set. rast n=1 
FFEA 57 FFFFFFFC 8F 54 Fi 0781 1198 3$ ACBL R4, #-4, R7, 2$ ; nn = until C 0) chowsh C(3) are in 
078B 110 ; though 
0788 1104 
0788 1105 ; 
078B 1106 ; Generate the low 128 bits of the product ity L. This product is 
0788 1107 ; equivalent to multiplying K times C'' modulo 8. The result of the | 
0788 1108 ; 13S it is in 14/17 with bits 31:29 of T4 the octant bits, and the remaining 
0788 1109: 35 bits the faction bits of the product. The last 53 fraction bits (bits 
078B 1110 ; 0 of T6 and 31:0 of T5) are non-valid fraction bits that will be used | 
078B 1111 ; ered if more fraction bits need to be generated. 
0788 Wie 3 
07868 1113; 
0788 1114 ; Multiply the high order bits of K (RO) times C'' and store the result in | 
0788 1115 ; TO/T2. 
078B 1116 
6— 00 04 A6 50 7A OQO78B 1117 EMUL RO, 4(R4), #0, (SP) 3; TO/T1 = KHI*C(3) | 
04 AE 04 AE 08 AG 50 7A 0791 1118 EMUL RO, 8(R4), 4(SP), 4(SP) ; he = KHI*(C(2):C€(3)] 
64 OC AG 50 cS 0799 1119 MULL3 RO, 12(R45, (R4) 3 14 = meh: bits if KHI*C(1) 
08 AE 64 C0 079 1120 DDL =—s-_« (RA), - BC SPS : 10/12 "2 "KHI*c'* modulo 8 | 
O7A2 1158 ; Multiply the low order bits of K (R1) times C'" and store the result in 
O7A2 1123 ; 1T4/T8. 
O7A2 1124 | 
64 00 04 A4 51 7A BoA 1125 EMUL R1, 4(R4), #0, (R4) 3 14/T5 = KLO*C(3) 
04 AG 04 A4 08 AG 51 7A O7A 1126 EMUL R1, 8(R4), 4(R4), 40R4) = T4/T6 = KLO*LC(2):C(3)) 
08 AG 08 AG «(0C AS «O51 7A—«(0780 «(1127 EMUL = R1, 12(R4S, B(R4S, B(R4): 14/1? = KLO*CC(1) :€(2):€(3)) 
51 C4 0788 1128 MULL R1, (R5) ; T8 = KLO*C(0) | 
65 CO 07BB 11¢9 ADDL = (R5), 12(R4) ; 14/17 = KLO®C'' modulo 8 | 
fet W13) 3; Add KHI*C'* to KLO®C’' to get K*C''. Store the result in 14/T7. 
co Me} 1138 ADDL (SP), 4(R4) 5 
D8 O7C3 1134 ADWC = &( SPS, B(R4) : 
dD8 O7C8 1135 ADWC 8(SP), 12(R4) : 1T4/T7 = K*C'* modulo 8 
O7CD 11 § 
O7CD 1137 ; : 
07CD 1138 3 At this potas there may or may not be enough, vel t¢ bits in R3/R4 to generate 
O7CD 1139; Y. irst 12 fratction bits ore all 1 r O's, there a possibility of 
07€D 1140 : Loss /* significance when computing, ¥ lanacquentts., we must check for loss 
O7CD 1141 ; of significance before converting 4/17 to Y and I. ; 
Br Ng 
C1 QO7CD 1144 ADDL3 #*X200000, -4(R5), (RS) ; If the first 12 fraction bits are 1's 
D3 0706 1145 BITL  #*x3FC00060, (R5) : and the reduced arg = 1-f or the 
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37. 12 «#2070D 1146 BNEQ CONVERT i first 7 bit are 0 and the reduced 
7DF 114 3 359 = f, then (and only then) bits 
70F 1148 : 9:22 are 0 and significance will 
7DF 1149 ; be lost. 
7DF 1150 
7DF 132) ; 
7DF 11 ¢ ; More bits need to be generated to cover the loss of Stent treanee. There are 
BepF 1155 ; not enough registers fo hold all the potential entra bits, so that the bits 
7DF 1154 ; already generated must be put on the stack. 
OPE 1186 ° 
QOOO0O90B'EF 16 O7DF 1157 JSB GEN_MORE_BITS ; Generate 85 additional bits and add 
O7E5 1158 3; =o them to oxigtiog bits. Results are 
O7E5 1159 ; . stored in 13/T7 
54 04 C2 O7E5 1160 SUBL #4, RG ; Adjust R4 to reflect the addition of 
O7E8 1161 ; _ another longword of K*C*’ 
15 FC AS 1D EO O7E8 1166 BBS #29, -4(R5), 4$ : Check if loss of significance is due 
1433 1187 ; to leading ones or zeros 
Bred 8 ; Lost significance due to leading zeros 
65 10A4 15 00 EA Q7ED 1167 FFS #0, #21, 16(R4), (R5) ; If at least one bit is set. This 
12 O7F3 1168 BNEQ CONVERT ; _ means lost significance was minor. 
OC AG TFFFFFFF 8F D1 O7F5S 1169 CMPL M°X1FFFFFFF, 12(R4) ; If one of the three high bits is set, 
17 15 «(O7FD «61170 BLEQ CONVERT : lost significance was minor. 
00B2 31 O7FF 1171 BRW LEADING_ZEROS 
0802 1106 2 é 
0802 BS ; Lost significance due to leading ones 
65 10A4 15 00 EB 0802 1175 4$: FFC #0, #21, 16(R4), (R5) |; If at least one bit is clear. This 
C 12 0808 1176 BNEQ C ; _.means lost significance was minor. 
OC AS €0000000 8F D1 O80A 1177 CMPL #*XE0000000, 12(R4) ; If one of the three high bits is. 
02 1& 0812 1178 BGEQU CONVERT ; clear, lost significance was minor. 
38 3611 «+0814 «1179 BRB LEADING_ONES 
0816 1180 
0816 1181 
0816 1186 CONVERT: 
0816 118 ; , 
0816 1184 ; Isolate octant bits and convert fraction bits to a pair of D-format 
pale 110) 3 quantities YHI and YLO 
65 FC AS 03 10 EF 0816 1187 EXTZV #29, #3, -4(R5), (R5) ; T8 = octant bits 
FC AS _E0000000 8F CA OB1C 1188 IcL_ #*x€0006000, -4(R5) : Clear octant bits ; | 
oc ¢3 base 1132 SUBL3 #12, RS, R4 7 R4 point to low order bits of h 
00000951'EF 16 0898 1190 SB «CVT. TO DOUBLE : RO/RI = 3s 9#h Lo | 
O82E 1191 3 R6/r7 = 2°29eh_hi | 
56 0680 8F A2 08 1136 SUBW #*XE80, R6 3 R6/R7 = h_hi | 
02 14 08 119 BGTR 3$ ; Check for h hi = 
564 1194 CLRL R6 ; Restore h_hi to 0 
50 83 1195 3$ TSTW RO ; Check for h_lo = 
05 1 9 1138 BEQL 1$ 
50 0€80 BF Ag 119 SUBW #*xE80, RO ; RO/R1 = h_lo ; 
07 65 9 0840 1198 18: BLBC =—-« (R55), 28 ; Check for odd or even octant bits | 
| 
1301 ; Octant bits are odd. Reduced argument equals 1 - h. 
56 08 56 63 43 1202 SUBD3 R6, #1, R6 3 R6/R7 = Y = 1 = bhi 
| 


; Increment weighting factor. If 
: weighting factor is greater than 
; 1024 then no more bits need to be 


| 

; Compress representation | 
Ke ee | 
generated. 
| 

| 


The og age factor is greater than 1024. This means that the reduced 
argument is either not distinguishable from 1 or too small to be represented 
in F-format (i.e. underflow.) Zero is returned in R4 for the reduced 
argument to signal this occurance. Note that under these conditions the 
correct function value is one of the values 0, +/-1. The . 
correc coetee is determined by the calling program based on the octant bits 
returned in R1. 


CLRL R3 ; Reduced argument is zero 
EXTZV #29, #3, 8B(SP), R2 ; R2 = octant bits 
BRW RESTORE 
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MTHSDSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VMS Macro V04-00 Page 28 
Bt O07 REDUCE_LA GE aes 07:59:38 UATHRTL. SREINTHDSINCO.MAR: 1 . (8) 
50 50 72 B47 1208 MNEGD RO, RO ; RO/R1 = -h_Lo 
2 : ‘ ; Get octant bits 
52. 20. A 0 084A 1207 28: MOVL 32(SP), R2 ; R2 = octant bits | 
Ooms a 4 1 08 BRW GET_YH{_YLO 
51 1 
51 1210 
31 1511 ; | 
51 1 \ ; At this point it has been determined that there is a major loss of 
851 1215 ; significance and the processing begins a looping phase. Each iteration of | 
p82) 1214 ; the loop will generate additional extra bits of K*C’ until enough significant 
51 1215 ; bits to compute Y are available. Goring this time the nine longwords 
tt ' 1§ ; allocated on the stack will be used as follows: | 
Bae) : 18 : TO/T2 Temporary storage used when generating extra bits. | 
8 3 : $0 : T3/T7 = Contains all significant bits generated so far. | 
0851 1222 : T8 Contains a counter, W, indicating the appropriate exponent 
0851 1558 ; of the Last longword of frection bite used in convert ine 
0851 1224 ; to Y. 
0851 1225 
S 3} : § LEADING_ONES: 
0851 1228 ; If processing continues here it is known that the loss of significance is due | 
8 2} : $3 3 to a string of leading ones. 
65 QO001E80 8F DO 0851 153) MOVL #L_INT_WEIGHT, (R5) ; «a oup bias for last longword 
8 28 1536 : of the product K*C' 
OC AS =6FFEQOOOO 8F) D1 s« 0858 1334 LOOP_1: CMPL #*XFFEQOO00, 12(R4) ; Check for enough significant bits 
2E 1A 0860 1235 BGTRU CONVERT 1 ; Enough bits. Convert to floating. 
0000090B'"EF 16 0 oe 1236 JSB GEN_MORE BITS : T2/T7 contains K*C'' 
OC A4 FFFFFFFF 8F D1 086 1237 CMPL #-1, 12¢R4) 3; Check for all 1's 
1E 1A 0870 1238 BGTRU CONVERT_1 ; Not all 1"s. Enough precision bits 
0872 1239 $ to compute Y 
08 A4 04 AG 7D 087 1240 MOVQ 4(R4), 8(R4) 5 
64 FC AG 7D 0877 1241 MOVa = =4( R45, (RG) . * 
FFD3 65 1000 8F 4000 8F 3D 0 78 4g ACBW #W_MAX_WEIGHT, #W_TERM_WEIGHT, (R5), LOOP_1 
0 iste 
0 : $2 
0883 1529 
0 1248 
0 1249 
0. 1250 
8 1251 
1 26 
0 125 
0 1254 
1386 
i 
1389 


0.09 0909 Cd 0 Cd 09 Gd C9 Cd 09 C9 GD Cd Od Cd OS Gd Cd Cd Od C9 Cd OD OS C9 Cd Cd CO CD GD CD CD CD 2.00 C2 COCO 
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8 
8 
8 
: 
8 
8 
8 
8 
8 
5 
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MTHSDSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VMS Macro V04-00 Page 29 
ait REDUCE. LARGE aes 94:59:38 UATHRTL. SREIMTHDSINCO.MAR: 1 ° (33) 
90 1260 
90 1261 
54 04 «60 3 : 34 rogers | #4, RG R4 poi to low bi fh 
3 oints to low bits o 
00000951'EF 16 0893 1564 JSB -CVf_T0_DOUBLE : RO/RI = 2*ueh_ lo 
99 1265 3 R6O/R7 = 2*Weh hi 
56 =—-FE72 gf 56 $3 99 1 96 SUBD3 =R6, D_2_T0_32, R6 3 R6/R7 = 2*ue(T = h_hid) 
50 SF 126 MNEGD ; RO/R1 = -2“*Weh Lo 
S52 FC AS OO 1p EF BAg 1 68 EXTZV #29, #3, -4(R5), R2 ; Re = octant bits 
54 $ C AB 126 SUBL (R55, Rb : RO/R7 = 1 = hh 
0 65 QO8AB 1270 TSTW ; Check for h_lo = 0 
4A 1 BAD 1271 BEQL GET_YHI_YLO : 
50 65 C AF 1 es SUBL ° ; RO/R1 = = h_lo 
45 11 0882 127 BRB GET_YHI_YLO 
08B4 1274 
08B4 1275 
0884 1 a LEADING_ZEROS: 
0884 127 
0884 1278 ; If processing continues here it is known that the loss of significance is due 
0884 1279 ; to a string of leading zeros. Note that it is known that the loop for 
0884 1280 ; leading zeros will terminate before an underflow condition occurs so that the 
bape ! + ; loop does not include a test for underflow. 
65 QO001E80 8F DO 0884 1 Hy MOVL #L_INT_WEIGHT, (R5) ; T8 = exp bias for last longword 
ores ! Be 3 of the product K*C* 
OC AS) OOIFFFFF BF D1 tra 1286 LOOP_O: CMPL #*XOOIFFFFF, 12(R4) ; Check enough fraction bits 
B 19 O8C3 1287 BLSS CONV i Enough bits. Convert to floating 
00000908 ‘EF 16 O8CS 1288 JSB GEN_MORE_BITS ; T2/R7 contain Ke*C** 
OC AS DS O08CB 1289 1$ TSTL 12(R4) ; Check for all O's , 
10 12 O8CE 1290 BNEQ CONVERT_O ; Not all O's. Enough precision bits. 
08 AS «6004 AG) «O7D SCOBD0 Ss 1291 MOVQ 4(R4), BC(R4) ; Compress representation 
64 FC AG 7D 0805 1 3¢ MOVa = =4(R4S, (R4) : | oof Kec’ | 
65 1000 8F AO Q8D9 129 ADOW #W_TERM_WEIGHT, (R5) 3: Increment weighting factor. 
DB «611 «#+O8DE 1294 8 LOOP 
O8E0 1295 
O8E 1596 CONVERT_O: : 
54 04 C0 pee 1297 ADD #4, RG : R4 points to low bits of h 
00000951"EF 16 O8E3 1298 CvT_TO_DOUBLE : RO/R1 = 2*Wth_Llo | 
O8E9 944 3 R6/R7 = 2*Weh_hi 
52 FC AS O03 1D EF O8E9 1 9 EXTZV #29, #3, -4(R5), R2 3 R2 =_octant bits | 
56 65 C2 OBEF 1301 SUBL (R55, Rb > R6O/R7 = h_hi 
50 65 O8F2 1 8 TSTW 0 ; Check for h_lo = 0 
t 1 pore 1 8 BEQL GET_YHI_YLO 3 
50 6 C th} : be SUBL (R53, RO ; RO/R1 = h_lo 
bare 1 8 GET_YHI_YLO: 
54 56 7D O8F9 1 vQ R6, R4 ; R4/RS = high bits of Y 
57 D& O8FC 1308 CLRL R7 ; R6/R7 = high 24 bits Y = YHI 
54 56 6 are 1309 SUBD R6, R4 : 
> BS ¢€ + 1310 ADDD RO, R4 3; R4/R5 = YLO 
904 1311 RESTORE: ; 
24 a 8 904 1 \¢ BISPSW 36(SP) ; Restore IV bit and exit 
C 907 131 ADDL #40, SP ; Remove mask and temporary storage 
0 90A 1314 RSB 
4 1315 
9 1316 
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MTHSOSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VMS Macro V04-00 P 
BOOP REDUCE. LARGE 66-SEP-1984 91:39:38 EMTHRTL. SREINTMDSINCO.MAR: 1 si 33) 
9 1317 GEN_MORE_BITS: 
Be 
908 1320 : This subroutine generates 85 extra fraction bits and puts them to the 
90B 1321 ; existing bits. NOTE: This routine is always entered via a JSB instruction. 
90B 1 ; 3 tet peace SP points to the first longword BEFORE 10, rather than T 
90B 1 3; itself. 
308 1358 
52 04 C2 0908 1 § SUBL #4, R2 ; Adjust polncer to get next quadword 
90E 1 ; rom MTH$AL_4_OV_PI 
56 62 33 79 395 1 8 ASHQ R3, (R2), R6 3: R7 = C(n) 
1 18 Bais : 3 BGEQ 1$ ; Branch if high bit is clear 
34h : 1 ; Logic to process unsigned values greater than 2°31 - 1 
04 AE 00 57 51 7A b914 1 é EMUL R1, R7, #0, 4(SP) s 
Ag 51 CO O91A 1334 ADDL R1. 8(§P) > TO/T1 = KLO®C(n) 
08 AE OB AE 5 50 A 34h; 1335 EMUL RO, R7, 8(SP), 8(SP) : 
OC AE 50 CO 09 1 § ADDL RO, 12(SP) : TO/T2 = K*C(n) 
OD 11 a : 38 RE 2$ 
89 : ' ib ; Logic to process unsigned values less than 2°31 
04 AE 00 57 51 7A 0928 1341 158: EMUL R1, R7, #0, 4(SP) ; TO/T1 = KLO*C(n) 
08 AE 08 AE 57 50 7A baa : g EMUL RO, R7, B(SP), 8B(SP) 3 TO/T2 = K*C(n) 
0338 : rf: 3; Add new bits to old 
64 08 AE CO 0938 1346 2$: ADDL B(SP), (R4) ; 
04 AG OC AE 08 093C 1347 ADWC =: 12 (SPS, 4(R4) : 
08 1E€ 9941 1348 BCC 3$ 3; Check for carry from previous add 
08 AS D6 0943 1349 INCL aa ; Propagate carry 
05 i€ base 1350 BCC $ ; Check for carry from previous add 
Oc AS 6 Soe 1351 INCL 12(R4) ; Propagate carry 
FC AG 04 AE DO 0948 1 26 3$ MOVL 4(SP), -4(R4) ; Move new low order bits to end of 
0950 135 : of old low order bits 
05 0950 1354 RSB 3 
348 1355 
951 1 2$ 
0951 135 
$98] 1399 
0381 1360 CVT_TO_DOUBLE: 
38} 1362 
$381 1 8 : This routine converts an array of three Longword pointed to by R4 to a pair 
951 1364 ; of D-format values. The results are returned in RO/R1 (low 48 bits) and 
2) : 9 ; R6/R7 (high 48 bits). ; 
50 Be 3 yy ' $ gyiLe | eee RO 3; RO/R1 = Low 32 bits of h 
07 «14 956 1 6 BGTR $ : Adjust for signed 
if D6 0958 1370 INCL (R4) 3 conversion error 
3 1€ QO9SA 1571 BCC 1$ ; If necessary, 
of AG bg 95C 1 ye INCL 4(R4) 3 __ propagate carry. 
50 1000 8F A 95F 1373 1$ SUBW #W_TERM_WEIGHT, RO ; RO/R1 = (low 32 bits of h)/2*32 


1 14 
MTHSDSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-19 1: AX/VMS v04-00 Page 31 
B00? REDUCE “LARG E g- Seen one 9 33: $ YRtARte. oATH HDSINCO.MAR; 1 . Fd) | 
FC AG 64 reer oone 8F C 4 1374 2%: BICL3 #*XFFFFOOOO, (R4), -4(R4) 
FC AG D 1375 SUBL.  =4(R4), (ROS : 
2 FC AG 1 16 CVTLD -4(R4), R2 ; 
50 g af ! 4 oper » RO ; RO/R1 = (low 48 bits of h)/2*32 
50 1900 F A2 O97A 1 8 SUBW #W_TERM_WEIGHT, RO ; RO/R1 = (low 48 bits of h)/2%64 
5 4 6 97F 1 9 3$: CVTILD (R&)+, R2 3; R2/R3 = next 16 bits of h 
9 1 982 1 BEQL 5$ 3 
¢ 14 0984 1 § BGTR 4$ ; aarees Te for signed cenyersten error. 
g + 9 § 1 INCL (R4) $ at no gerry is poss {dhe 
52 1000 8F rl 9 1384 4$ SUBW #W_TERM_WEIGHT, R2 3 R2/R3 = "ieeen 16 b of h)/2*32 
56 64 «4966 «098D 1385 5$ CVTLD (RZ), RB : Ré/R? = high 32 bits of h 
By 9 1 $ BGEQ 3 a =P for signed conversion 
56 FO7A CF 6 992 1 : ADDD D_2_T0_32, R6 : 
56 2 60 0997 1388 6$ ADDD R2, R6 : RO/R? = = (high 48 bits of h)/2*32 
a: et 
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Es -1984 11: MTHRTL.S SMTHDSINCO. MAR; 1 


-SBTTL REDUCE _DEGREES 


: This routine assumes that the seatyte value of the argument is in RO/R1. 
: The reduction process is performed in two gigges. The first stage 


nt Sine, Cosine ond OXnena y 2 tae + 9 1: +30: 3 ales Be o v04-00 Page | 


the reduction reduces the argument modulo value less that 2°55, 
; and the second stage reduces the argument modulo és to a value less than 45. 


: Constants used in this reduction: 


POWER HOD, 360.0: , ; Powers of 2 modulo 360 for t1 = 0 
“WORD 6-16, 3 - 64, 12 
«WORD 256, 152, 304, 248 
POWER _MOD 360. 1: ; Powers of 2 modulo 360 for t1 <> 0 
.QORD™ 136, 2872, 4, 8 
“YORD 16, > Ste 128 
«WORD 256, 152, 304 248 
REDUCE_DEGREES: 
” CMPW #*x5C00, RO : Compare ix! with 2°55 
BGTR LAST_S TEP ; Branch to special logic for med arg 


It is assumed here that the argument is greater than 2°55. 


The argument is reduced as follows: 

ef x = 2*tef, where t ; 56 and 1/2 =< f < 1. And let J = 2°56*f = 

“ Ory1 + J2 and K = 2*(t-56). Since 236 : 64 modulo 560, we have that 
J_= 64*j1 + J2 modulo 36 5. Now let .* =t - 56= ats tl + t2. Note that 
(2*12)*2 = (2°9)#(2°15) = (2*9) #(2°3) £43 modulo 360. Hence, if ti is 
not zero, K = 2*t* = 2°(12*t1+4t2) = (2"1 )* (2*t2) = 136% 2*t3 agdylo 700. 
For tl = 0K = 2*t2. Consequently, define Kk‘ congruens to 2*t2 if tl = 0 
and congruent to 156*2*t2 otherwise, where 0 =< < 360. Then x‘ = 
K'#(64*J1 + J2) is congruent to s modulo 360 cae sg ¢€ oh 


aa aaa a ed a tt wm 


PS ee ee ee ee ee ee ee ee ara 
ONAN E WO ODNOANE WIN OOO NOU EWN 0 ODNAUE WN OOD NOAANEWN OOO NOU EU 


BBE BE EE EMI ANIA IPP NPINPININININID) 2 9 S| 2 QOOOOOOOCOOCOWOOOOoOowowo 


MOVL Re ; R2 = high longword of X 

BICL Pare F80, RO ; Clear exp bits of X 

B1SW tds ae RO ; RO/R1 = J 

SUBL RO, R2 ; R2 = t'#2*7 

mova RO, R3 3; R3/R4 = J 

BICL #°XFFERSFFF, R1 : RO/R1 = J1*2*30 

SUBD ° : R3/R4 = J2 

SUBW #*xC00, RO 3 RO/R1 = 64*J1 

ADDD R3, RO ; RO/R1 = 64/1 + J2 = J modulo 45 
ROTL #-7, : - Re 3; R2 = t' 

DIvw3 ne. k R ; R35 = tl 

MULW ; RS = tal 

SUBW ; R2=t 

TSTW = RB. ; Check for tl = 0 and choose K’ 
BNEQ 1$ 3; =©accordingly 
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MTHSDSINCOS ; Floating Point Sine, Cosine wt Sincos a EP=1984 01: AX/VMS v04-00 Page 33 
Brb0F REDUCE. DEGREES Breer e8s FT iSki8s LATA NTS eects oer OO mans1 29 B, 
52 BB AFK2 GD ane 1449 CVTWO — POWER_MOD.360_0(R2I, R2 ; R2/R3 = K* | 
52 45 AFG 6D QAI 1681 1$: CvTwDd + a MOD_360_1CR2], R2 : R2/R3 = K' 
5 Al 14 § 2$: MULD R2, RO 3; RO/R1 = X" (mod 45) 0 =< RO < 2°55 
Ag Ue | 
rt 142? LAST STEP: 
aie 1437 : Argument reduction scheme for arguments with absolute value less than 2°55 
AiB 1238 : The reduced argunent Y is computed as follows: 
AiB 1299 ; Let I = int(x/45) 
AIB 1461 ; if s even 
OA1B 1096 3 then Y = X = 45* 
QOAIB 14635; else Y = (141)#4 
OAIB 1468 
50 5200 of * OAs 1298 oy onpenes, oe - ; Compare 2°36 with ‘xX! 
56 50 F622 CF 65 OA ; 1468 MULD3 D_T_O OV wis RO, R6 : R6/R7 = {x!/45 
bA58 1490 
- f \o0) : Turn off IV to avoid an exception in EMODD 
52 DC OA28 1298 , MOVPSL R2 ; Move PSL to R2 
52 FFFFFFDF 8F CA - A 1674 BICL #*C<PSLSM_IV>, R2 ; Save current IV bit 
20 «Bg ny 132? BICPSW #PSLSM_IV ; Turn off integer overflow trap 
54 53 56 00 O08 7% Saas 1477 EMODD #1, #0, R6, R3, R4 ; a low 32 integer bits of ixi/45 
OARS 1428 : R4/RS = fractional part of 'x:/45 
52 B68 one rt + BISPSW R2 3; Restore IV bit 
6,, 33 6¢ pase 1482 sees Ree even ; R6/R7 = Integer part of :xi/45 = I 
56 08 60 Saas 1484 ADDD2 #1. R6 : R6O/R7 = 14+ 1 
12 11 OA44 1485 BRB ODD 
ae 1s 
A46 1488 NO_OVERFLOW: 
54 53 50 0B FSFE CF 74 we 1932 EMODD 0_1_0V_45, #X_1_0V_45, RO, “a, 3 etenie ar ah ekiee 
1453 €9 OASE 1491 BLBC Ree CvT : Branch if oc that its are even 
56 33 9} C1 vel 1636 ADOL3 R3, R6 3; R6 =] + +, 
6 6 6E OASS 149 CVTLD 48 R6 3 R6/R7 = 
+ sok Se ME: EY): ena 3 A a t RUA? éseclet) 
53 FB BF BA A60 1496 BICB. = s #*KFB, R3 : Save only last three octant bits 
> Aes 1298 ~— 
6 53 6& OA65 1499 CVT: CVTILD =R3, R6 3 R6/R7 = I | 
76 EBC ot 8S Bago 1301 Sven: ae BAG, Re 4+ + Pees 
53 ~F8 8F bh A70 13 § BICB #°XFB, R3 > Save only last three octant bits | 
5 Oa74 15 RSB | 
A75 1504 
A75 1505 | 


; 
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MTHSDSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VMS Macro V04-00 Pa 4 
ah 7 REDUCE. DEGREES or see Oke 91:39:38 EMTHRTL. SREIATHDSINCO.MAR: 1 ” (33) 
A75 1507 
Ae } os -SBTTL RADIAN_POLYNOMIALS ; Polynomials for arguments in radians 
A75 1510 
A7 1511 
A> 1518 | 
ns 131? ; Polynomial evaluation for DCOS(Y) for Y in radians | 
A7 1313 
A?5 1 16 P_cos_r: 
53 29 8000 8F AB OA75 151 BICWw3 #*x8000, Rg. R3 3 
5 4000 8F BI 4 1318 CMPW #*x4000, R 3; Compare 1/2 with !YHI! 
1 14 QOA8O 153 BGTR LEQL_HALF ; Sufficent overhang is available 
OA8 1520 NEEDS_DOUBLE: 
7E 54 7D me 1521 MOVQ R4, -(SP) ; Save YLO 
7E 546 56 61 A 15 § ADDD3 = =R6, R4, -(SP) 3; Save Y 
50 6e gE $3 Ago 15 MULD3 (SP), (SP), RO : RO/R1 = Y 
F675 CF 0 5 5 gaat 1524 POLY RO #COSLENR2-1, COSTBR2; RO/R1 = Q(Y*2) 
54 56 BE 61 0A93 1525 ADDD3 (SP)+, R6, R4 : R4/RS = Y + YHI 
54 BE 64 OA97 1526 MULD (SP)+, R& 3 R4/RS = YLO#(Y + YHI) = A2 
05 13 OA9A 1527 BEQL : Check for A2 = 0 
54 0080 8F As OA9C 1528 SUBW #*x80, R4 : R4/R5 = A2/2 
50 54 62 OAA1 1529 1$ SUBD =—s_- RG, -RO + RO/R1 = Q(Y*2) = A2/2 
56 56 64 OAAS& 1530 MULD R6, R6 : R6/R7 = YHI* 
56 0080 8F A2 OQAA7 1531 SUBW #*x80, R6 3; R6/R7 = YHI*2/2 
56 08 62 OAAC 15 : SUBD = #1, _ RO : R6/R7 = -(1 = YHI*2/2) 
50 56 6 OAAF 15 SUBD R6, RO ; RO/R1 = DCOS(Y) 
05 OAB2 1534 RSB 
OABS 1535 
OAB 1336 LEQL_HALF : 
56 54 60 OAB3 1 $ ADDD R4, R6 ; R6/R7 = 
56 56 64 OAB6 1538 MULD R6, R6 3; R6/R7 = Y*2 
F609 CF 07 56 75 OAB9 1539 POLYD R6, #COSLENR1-1, COSTBR1; RO/R1 = DCOS(Y) 
05 OABF 1540 RSB 
QACO 1541 
OACO 1346 
OACO 1543 ; 
OACO 1544 ; Polynomial evaluation for -DCOS(Y) 
OACO 1545 ; 
OAC 1366 
AC 1547 N_COS_R: 
53 56 8000 8F AB OACO 1548 BICWS #*x8000, R6, R3 3 
53 4000 8F B61 OAC6 1549 CMPW #*x4000, R3 ; Compare 1/2 with ‘YHI: 
2 14 OACB 1230 BGTR ; Sufficent overhang is available 
7E 4 7D OACD 1551 MOVQ R4, -(SP) ; Save YLO | 
7E 54 56 «61 ADO 1336 ADDD R6, R4, =-(SP) 3 Save Y 
50 «6E gf $2 ADS 155 MULD3 (SP), (SP), RO : RO/R1 = Y 
F62A CF 2 5 5 AD8 1554 POLY RO ACOSLENR2-1, CCOSTBR2; RO/R1 = Q(Y*%2) 
54 56 BE 61 OADE 1555 ADDD3 (SP)+, R6, R4 > R4/RS = Y + YHI 
54 E 64 AE 1338 MULD (SP)+, R4 3: R4/RS = YLO®(Y + YHI) = A2 
13° OAES 155 BEQL 1$ 3: Check for A2 = 
54 080 8F as AE7 1338 SUBW #*x80, R4 3 R4/R5 = A2/ 
0 54 62 OAEC 1559 18: SUBD R4, RO : RO/R1 = Q(Y*2) = A2/2 
6 6 64 OAEF 1560 MULD R6, R6 ; R6/R7 = YHI* 
56 080 8F A AF 1561 SUBW #*xX80, R6 : R6/R7 = YHI*2/2 
4 8 62 OAF 136¢ SUBD. = #1, - Rb > RO/R7 = -(1 = YHI*2/2) 
50 6 6 AFA 156 SUBD3 RO, R6, RO ; RO/R1 = -DCOS(Y) 


————— nn 
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MTHSDSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01: AX/VMS 0 V04-00 Page 35 | 
ity, RADIAN. POLYNOMIALS ; Polynomials for arg 6-SEP-19 1986 ti 38: Hy EMTHRTL. SRE IMTHDSINCO. MAR; 1 ° (33) | 
0S OAFE 1564 RSB | 
AFF 1565 
56 54 60 QAFE 1566 2s: ADDD RG, - RG : R6/R7 = Y | 
36 6 64 : 156 MULD 6, R6 + R6/R7 = Y*2 
F SBD cf 0 36 5 1368 POLYD R6, #COSLENR1-1, COSTBR1: RO/R1 = DCOS(Y) 
0 8000 8 ac OB0B 1 : XORW  #*8000, RO + RO/R1 = =DCOS(Y) 
05 0810 1570 RSB 
111571 | 
B11 1326 ; 
B11 1573 ; Polynomial evaluation for -DSIN(Y) 
11. 1574 : 
Bi 1323 
161 i: N_SIN_R 
54 8000 BF AC O0B11 157 "YORW -#*XB000, Ré ; 
56 8000 BF AC 0816 1578 XORW  #*X8000, R6 t RG/R7 = -Y 
BiB 1579 
0818 1580; 
0B1B 1581 ; Polynomial evaluation for DSIN(Y) 
0818 1386 3 
0B1B 15 
0B1B 1384 P_SIN_R 
7E 54 7D OB1B 1585 MOVQ = R&, : Save YLO 
7E 56 54 61 OBIE 1586 ADDD3 RG, RE, - : Save 
54 6 6€ 65 082 1587 MULD3 (SP), (SP), R4 : RG = Y* 
FOIC CF «3600? 54 SOS (B26 1588 POLYD R4, #SINLENR-1, SINTBR : RO/R1 = P(Y*2) 
0 BE 64 OBec 1589 MULD (SP)+, RO + RO/R1 = Y*P(Y*2) 
0 8 60 OB¢F 1590 ADDD (SP)+, RO + RO/R1 = YLO + Y*P(y*2) 
0 56 60 OB32 1591 ADDD ~—R6, RO + RO/R1 = Y + Y*P(Y*2) = DSIN(Y) 
05 0B35 1398 
0B36 159 
0B36 1594 
0B36 1595 
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Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VME Macro V04-00 Pa 
EVCLE POLYNOMIALS : Polynomisle for orgy GoSEPH1ORG $4555585 FANANTE  eeeSerMocioco.mans1 2% (39) 
: 6 1338 -SBTTL CYCLE POLYNOMIALS ; Polynomials for arguments in cycles 
o%6 1398 | 
1600 
B36 1601 ; 
: 6 1 § :; Polynomial evaluation for DCOS(Y) for Y in cycles | 
B 4 1604 
6 1605 P_COS_C: 
56 F4GEE CF 71 B36 1906 CMPD D.2_0V_PI, R6 3; Compare 2/pi with {YHI! 
fo. B 16 BGEQ 23 : Sufficent overhang is available 
7 4 7D p 16 : MOVQ R4, =(SP) : Save YLO 
7E 5 6 61 BS 16 ADDD R6, R4, -(SP) 3; Save Y 
50 6e gf 65 0B44 1610 MULD3 (SP), (SP), RO : RO/R1 = Y 
F67A CF 2 75 0B48 1611 POLYD RO ACOSLENC2-1, COSTBC2; RO/R1 = Q(Y*%2) 
54 56 «BE 61 (0B4 1618 ADDD3 (SP)+, R6, RS : RG/RS = Y + YHI | 
54 a 64 O0BS52 161 MULD (SP)+, R4 3 R4/RS = YLO*(Y + YHI) = A2 
0 13° 0B5 1614 BEQL 3; Check for A2 = 0 
54 0100 8F as 0B57 1615 SUBW #*x100, R4 3 R4/RS = A2/4 
50 54 62 OBSC 1616 1$ SUBD R4, RO s RO/R1 = Q¢(Y¥%2) = A2/4 
56 56 64 OBSF 1617 MULD R6, R6 : R6/R7 = YHI*2 
56 0100 8F A Babs 1618 SUBW #*xX100, R6 3s R6/R7 = YHI*2/4 | 
56 8 62 086 1619 SUBD #1, R6 3 R6/R7 = =-(1 = YHI*%2/4) | 
50 56 62 OB6A 1620 SUBD R6, ; RO/R1 = DCOS( 
05 0B6D 1621 RSB | 
OB6E 16 ; 
56 54 60 OB6E 1625 2$ ADDD R4, R6 3 R6/R7 = Y 
56 56 64 0B71 1624 MULD R6, R6 3: R6/R7 = Y*2 
F6OE CF 07 56 75 0B74 1625 POLYD R6, #COSLENC1-1, COSTBC1; RO/R1 = DCOS(Y) =- 1 
50 08 60 OB7A 1626 ADDD #1, RO ; RO/R1 = DCOS(Y) 
05 087D 1627 RSB 
SB7e 1628 
7E 1668 
OB7E 1630 ; 
0B7E 193) ; Polynomial evaluation for -DCOS(Y) 
OB7E 16 § : 
OB7E 163 
OB7E 1634 N_COS_C 
56 FAG CF 71 pare 1635 CMPD D.2_OV_PI, R6 ; Compare 2/pi with !YHI! 
32.18 83 1636 BGEQ 2$ : Sufficent overhang is available 
7E 54 7D 0B85 1937 MOVQ R4, =-(SP) 3; Save YLO 
ct SS %- pase 1638 ADDD3 =R6, R4, -(SP) : Save Y 
50 6E of 65 OB8C 1639 MULD3 (SP), (SP), RO : RO/R1 = 
F632 CF 2 5 75 90 1640 POLY RO ACOSLENC2-1, COSTBC2; RO/R1 = Q(Y*2) 
4 56 B8E€ 61 0896 1641 ADDD3 (SP)+, R6, R4 > RG/RS = Y + YHI 
54 ge 64 OB9A 164¢ MULD (SP)+, R4 ; R4/R5 = YLO®(Y + YHI) = A2 
13 OB9D 164 BEQL 3; Check for A2 = 0 
54 0100 8F as OF 1644 SUBW #*xX100, R4 3 R4/RS = A2/4 | 
50 4 6 A4 1645 1$ SUBD R4, RO ; RO/R1 = Q¢Y*%2) = A2/4 
56 56 64 A7 1948 MULD R6, R6 3; R6/R7 = YHI* 
56 100 8F A AA 164 SUBW #*xX100, R6 : R6/R7 = YHI*2/4 
6 08 6 AF 1968 SUBD ‘ 3 R6/R7 = -(1 = YHI*2/4) 
50 56 50 6 BB2 164 SUBD3 RO, R6, RO 3; RO/R1 = -DCOS(Y) 
> Ope eat ans | 
56 54 60 0BB7 1926 2$ ADDD R4, R6 3; R6/R7 = Y 
56 56 64 OBBA 165 MULD R6, R6 3; R6/R7 = Y*2 
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MTHSDSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VMS Macro V04-00 a 7 
a EYCLE POLYNOMIALS ; Polynomials for argu on 3eE=18 4 91:99:38 EMTHRIL. SRCINTHDSINCO.MAR: 1 ses (3) 
F5C5 CF 7 7 BB 1654 POLY R6, #COSLENC1=-1, COSTBC1; RO/R1 = DCOS(Y) = 1 
50 00000900 0000C080 7 28 63 BC3 1655 SUBD Re. #-1, RO 3; RO/R1 = -DCOS(Y) 
BCF 1656 RSB 
BDO 1657 
BD 1928 3 
BDO 1659 ; Polynomial evaluation for -SIN(Y) 
BDO 1660 ; 
BDO 1661 
BD 1006 N_SIN_C: 
54 8000 8F AC OBDO 166 XORW #*x8000, R4 3 
56 8000 8F AC OBD 1664 XORW #*x8000, R6 3 RG/R7 = - Y 
BDA 1665 
BDA 1666 ; 
OBDA 1667 ; Polynomial evaluation for DSIN(Y) 
OBDA 1668 ; 
OBDA 1669 | 
OBDA 1670 P_SIN_C 
7E 54 7D OBDA 1671 MOVQ R4, -(SP ; Save YLO | 
7E 56 54 61 OBDD 1676 ADDD R4, R6, -(SP 3; Save 
54 6f 6 65 OBE1 167 MULD3 (SP), (SP), R& t RG = Y*2 
F61D CF 0 54 75 OBES 1674 POLYD R4, ASINLENC-1, SINTBC ; RO/R1 = P(Y*%2) 
50 BE 64 OBEB 1675 MULD  (SP)+, RO + RO/R1 = Y*#P(Y*2) 
50 6E 60 OBEE 1676 ADDD (SP), RO : RO/R1 = YeP(¥*2) + YLO 
6E 0100 8F A2 OBF1 1677 SUBW 3s #*x 160, (SP) : (SP) = YLO/4 
50 BE 62 OBF6 1678 SUBD (SP)+, RO : RO/R1 = YeP(Y¥*2) + 3/4*YLO 
54 56 7D OBF9 1679 MOVQ R6, R 3 R4&/R5 = YHI 
54 0100 8F A2 OBFC 1680 SUBW #*xX100, R4 ; R4/RS = YHI/4 
56 54 62 0CO1 1681 SUBD R4, R6 3: R6/R7 = 3/4*YHI 
50 56 60 OCO0& 1682 ADDD R6, RO ; RO/R1 = DSIN(Y) 
05 OCO7 1683 RSB 
0CO8 1684 


BT IRpSENCOS 
2-00 


56 F440 4 

54 6 

FOOF CF of & 
54 56 00070000 gf 
52. 56054 
5654 

56 3 

0 

56 0680 8F 

50 56 

5454 

54 0680 8F 

54 (08 

5054 

56 56 

07 

FSF7 CF 607) «6556 
50. (08 

56 FSFE CF 

38 

54 56 856 

F621 CF 07 = 54 
54 56 00070000 of 
a6 6 32 

560 54 

56 ¢ 

0 

56 0680 8F 

a 26 

4 4 

54 0680 8F 
5408 

50 54 50 

56 56 

9¢ 

FSA8 CF 07 6 
50 8000 8F 
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pine. Cosine and Sincos 16-SEP-19 


6-S SEp= 1984 i 39: 33 ic 


MTHRTL.SR 


-SBTTL DEGREE_POLYNOMIALS 


90_OV_PI, R6 


R6 
RCSSDINI cospret 
0000, 


#*xX680, R6 
Be. ey 


ari80, RG 
#1, 
R4, 0 


R6, R6 
R6, #COSDLN2, COSDTB2 


#1, RO 


D.90_OV_PI, R6 
2$ 


R6, R6, 

ae AcosbuNt COSDTB1 
: x70000, R4 

R4, R6, R2 

R4, R6 

R2, R6 


#*X680, R6 
R6, R R? 


R4, 
#ri80, R4 
#1, 
RO, Ra, RO 
R6, R6 
33° 
R6, #COSDLN2, COSDTB2 
#*x8000, RO 


-1, RO 


X/VMS Macro V04-00 


€JMTHDSINCO.MAR; 1 (24) 


Compare 90/pi with Y 
ReuR se ore ision isn't needed 


; RO/R1 acy*2) 

; R4/R5 = YHI 

; a YLO 

3; R6/ /RT = YHI 

3 R6/R7 = YLOsty + YHI) = A2 
3; Check for A2 =_0 

3 R6/R7 = “a 

3; RO/R1 = Q(Y*%2) = A2/2*13 
; R4/RS = YHI* 

3 R4/RS = YHI*2/2°13 

3 R4/R5 = =-(1 = YHI*2/2*13) 
; RO/R1 = DCOS(Y) 

3 R6/R7 = Y*2 

3; Check for Y = 0 

; RO/R1 = Q(Y*%2) 

3; RO/R1 = DCOS(Y) 


Compare 90/pi with Y 
pous : ores ‘ston isn't needed 


RO/R1 = QcY*2) 


ey YHI 
fF 
RG/R? = YHI 
R6/R7 = WLosty + YHI) = A2 
Check for A2 = 
R6/R7 = A2/2*13 
RO/R1 = A ) = A2/2*13 
R4/R5 = = 
R4/R5 = YHI"5/ 2°13 
R4/R5 = -(1 = Mitel 
RO/R1 = -DCOS(Y 
R6/R7 = Y*2 
Check for Y = 0 
RO/R1 = DCOSD(Y) 
RO/R1 = -DCOSD(Y) 
RO/R1 = DCOS(Y) 
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MTHSDSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:20: AX/VMS Macro V04-00 Page 39 
BrHeD BEGREE. POLYNOMIALS ae a 07:59:38 EMTHRIL. SREIMTMDSINCO.MAR: 1 . (a) | 
05 OCcB1 1743 RSB 
CB2 1744 
CB2 1745 N_SIN_D: 
56 56 72 OCB 1746 MNEGD R6, R6 3 R6/R7 = -¥ 
CBS 1747 P_SIN_D: 
soe eee Be cB5 1f03 MULD3 R6, R6, RO ; RO/R1 = Y*2 
11 13 OCB9 174 BEQL RETURN 
F607 CF 97 20 7 CBB 1739 POLYD RO, #SINDLN, SINDTB 3 RO/R1 = PC(Y*2) 
6 64 OCCI a 1 MULD . ; RO/R1 = Y*P(y¥*2) 
56 2 00 ef A ccé «(1 § SUBW #*X300, R6 s R6/R7 = Y/2°6 
6 g cc9 (17 ADDD R6, R s RO/R1 = DSINC(Y) 
33 Vf ¢ RETURN: RSB 
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2-007 


50 (08 


50 00000000 0000C080 8F 


a we 
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SQ 
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= 


6-SEP=19 


~SBTTL DEGENERATE SOLUTIONS 


MOVD #1, RO 
RSB 


RSB 


; Underflow; if user has FU set, signal 


MOVPSL R2 
CALLS aye G*MTHSSJACKET_TST 


BLBC 
MOVZWL SF. SAVE_PSW(FP), R2 
CLRL 


BBC #6, R2, 208 
PUSHL 
MOVZBL #MTHSK_FLOUNDMAT, -(SP) 


CALLS #2, G*MTHS$SIGNAL 
RSB 


~ END 


4 11: 


:52:38 EMTANTL. SREIATHDSINCO.MAR: 1 


; Answer is 1 


; Answer is -1 


error. Always return 0.0 
; ; R2 


; branch if user did 
et user PSL saved b 


x 
; has wee enabled floatin 
yes, return PC from spec 
: tres code for hardware floa 


; signal (condition, P 
; return 


user's or jacket routine’s PSL 
RO = TRUE if JS eer J jacket routine 


; RO = result. CIBSSIGNAL wi will save in 
; ox ~ “paetenian so any handler can 


underflow convert to MINS. Mt UNDMAT 
(32-bit VAX-11 exgget ies Code) 
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Symbol table 
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Oo 
wv 
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MTHSDSINCOS 
Symbol table 


Dlg 
SMALL =£080 
MA 


L.CcOos 


JUST 
uoRAX WEIGHT 
: TTERA_WE IGHT 


x ~1.0V_45 


Initialization 


ng processing 
panes, table sort 
Symbol table output 
Psect synopsis output 


Cross-reference output 
Assembler run totals 


The ya sock ing ¢ 


+ bong 
there were 


Macro Library name 


G 15 
; Floating Point Sine, Cosine and Sincos 16-SEP-19 


Allocation 
00000000 < 
00000000 <( 
00000001 


Page faults CPU 


is GR:aR-aB-og 9808: 
206 0:00: 3 00:00: 
0 3 B08°§ 00:00: 
320 9:00: °90 90:00: 
SBS Bee 
2 80:00; 8°08 00:00: 
69 00:00:11.02 00:00: 


t Limit was 1650 pages. 
6 pages) of virtual memory were used to buffer the intermediate cod 
0 pages of symbol table peese gl tgcated to hold 194 non-local and 57 lecat symbols. 
850 source 2? jpeg were read in Pass = 

15 pages of virtual memory were used to define 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 
131 GETS were required to define 5 macros. 


producin 


: 
| 


guosececceccce 


PSECT No. 


0¢ 0.) 


Attributes 

NOPIC USR 

NOPIC USR 
PIC  USR 


See OO EEE OOO ODM em + 


mOoONMronrnr 


object records in Pass 


macros. 


were eeeoeeeeeoeeeeeeeeseeeces} 


6-SEP-1984 


CON 
CON 
CON 


Be 91:82:38 


ABS LCL 
ABS LCL NOSH 
REL LCL 


AX/VMS M 
MTHRTL. SRE 


SHR EXE 


o V04-00 
JMTHDSINCO.MAR; 1 


Page 


NOSHR NOEXE NORD NOWRT NOVEC BYTE 
HR EXE RD WRT NOVEC BYTE 
RD NOWRT NOVEC LONG 


( 


42 
36) | 


v4 


TISEREIMESS, un statistics Fating Point Sine, Cosine and'stncos 1g-568>1364 1:88:38 ERTANTE. SREGRrNDSiNcO.man: 1 P29" (58) 


There were no errors, warnings or information messages. 
MACRO/ENABLE =SUPPRESS1ON/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDSINCO/OBJ=OBJ$:MTHDSINCO MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) MS 
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